5 APR 1949 

NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


» 

) 

TECHNICAL NOTE 
No. 1436 

» • 

STRESSES IN AND GENERAL INST ABILTT Y 
OF MONOCOQUE CYLINDERS WITH CUTOUTS 
VI - CALCULATION OF THE BUCKLING LOAD OF CYLINDERS 
WITH SIDE CUTOUT SUBJECTED TO PURE BENDING 
By N. J. Hoff, Bertram Klein, and Bruno A. Boley 
Polytechnic Institute of Brooklyn 



Washington 
March 1948 


— sen 

-' v A C A LiiikAKI' 

WElfOWAL 4i 
L/JSORA'ro*V 

frlftW, T- 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE NO. 143 6 


STRESSES IN AND GENERAL INSTABILITY 
OF MONOCOQUE CYLINDERS WITH CUTOUTS 
VI - CALCULATION OF THE BUCKLING LOAD OF' CYLINDERS 
WITH SIDE CUTOUT SUBJECTED TO PURE BENDING 
By N. J. Hoff, Bertram Klein, and Bruno A. Boley 


SUMMARY 


A strain-energy theory was developed for the calculation of the 
"buckling load in general instability of circular reinforced monocoque 
cylinders having a side cutout and subjected to pure "bending. The 
theory was applied to two series of specimens., each: containing three 
cylinders, which, were tested and reported in part TV in this program 
at the Polytechnic Institute of Brooklyn Aeronautical Laboratories. 
The average deviation between theoretical and experimental buckling 
load was 27.1 percent for the first series and 3^.4 percent for the 
second. 


INTRODUCTION 


In the present report, the last report of a series of six (see 
references 1 to 5) dealing with monocoque cylinders with cutouts, the 
calculation of the buckling load of a cylinder which has a. side cutout 
a n d fails in general instability when subjected to pure bending is 
undertaken. General instability is defined as the simultaneous buckling 
of the circumferential and longitudinal reinforcing elements of a mono- 
coque cylinder together with the sheet attached to them. As the cal- 
culations given herein follow closely those presented in reference 3 > 
it should be consulted for the development of some of the fundamental 
connections used in the present report . 

The first step in the strain-energy calculation was the assumption 
of a buckled shape. This was done after examination of the d eflecti on 
patterns observed in the tests described in reference 4. The stringer 
bordering the cutout on the compression side always showed the greatest 
deflections, and its distorted shape was very similar to a full sine 
wave . For the specimens having a symmetric cutout (reference l) the 
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wave length was only slightly greater than the length of the cutout, hut 
for the specimens with side cutout the length of the wave varied and was 
observed to he definitely greater. For this reason, in the present report 
the wave length L was considered a parameter, the value of which had to 
he determined from the requirement that the buckling load he a minimum. 
This procedure differs from that adopted in the calculation of the general 
instability of the cylinders with symmetric outout in which the wave 
length was assumed to he a known constant, namely the length of the 
cutout. 

In the circumferential direction the deflected shape at - buckling 
is represented by the first-seven terms of a Fourier expansion. The 
circumferential coordinate is measured from the edge of the cutout, and 
the length of the Interval in which the Fourier series Is defined Is 
considered as one of the parameters of the problem. The boundary con- 
ditions at the end of the interval determine four of the seven coeffi- 
cients of the series, and one of them is indeterminate as in all buckling 

problems. The remaining two coefficients, as well as the wave-length 

parameter, are calculated from the requirement that the buckling load be 
a minimum. 

The assumptions made regarding the buckling pattern in the circum- 
ferential direction are identical for the cylinders having a symmetric 
cutout and for those having a side cutout. For cylinders having a side 
cutout, however, additional considerations were necessary because the 
deflections extended Into the complete portions of the cylinder. It' 
was decided to assume that - the expressions developed for the buckled 
shape in the cutout portion were also valid for the complete portions 
of the cylinder J this assunption, in effect, means that the restraint 
due to the continuity of rings and sheet in the complete portions Is 
neglected. The Justification for this assumption Is the resulting com- 
parative simplicity of the calculations, as well as the observation that 
the deflections always were considerably smaller in the complete portions 
than in the cutout portion. 


The following strain-energy quantities were considered: radial 

and tangential bending, as well as torsion of the stringers j bending of 
the rings in their plane; and shear in the sheet. The extensional 
strain energy In the sheet was taken into account by adding an effective 
width of sheet to the stringers and the rings. In the calculation of- 
the external work it was assumed that the applied moment' caused a linear 
distribution of the strain in the cutout portion of the cylinder. The 
forces corresponding to these strains were applied to the stringers at 
the ends of the wave, and the sum of the products of these forces and 
the displacements of their points of application was taken as the 
external work . 

The buckling load was calculated from the requirement that the 
strain energy corresponding to the transition from the unbuckled into 
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the "buckled shape he equal to the work done hy the applied loads . The 
minimum value of the buckling load was found hy assuming the circum- 
ferential wave length to he equal to the length of some integral number 
of stringer fields, and the axial wave length, some integral number of 
ring fields. The values of the two undetermined Fourier coefficients 
were calculated to make the buckling load a minimum. This minimum value 
of the buckling load then was determined and compared with other minimum 
values obtained on the basis of different choices of circumferential and 
axial wave lengths. Between 6 and 10 combinations were investigated for 
each of the six cylinders - three of which are shown in figure 1 - in 
order to find the absolute minimum value of the buckling load. 

The investigation was conducted at the Polytechnic Institute of 
Brooklyn Aeronautical Laboratories under the sponsorship and with the 
financial assistance of the National Advisory Committee for Aeronautics. 
For his substantial share in the numerical work the authors are indebted 
to Mr . John G. Pulos. 


SYMBOLS 


Fourier coefficients 

cross-sectional area of stringer plus its effective width 
of sheet 

Fourier coefficients 

geometric factor in torsional rigidity GC 
width of panel measured along circumference 
Young's modulus 
shear modulus 

shear modulus of sheet covering at zero compressive 
load 

G 0 ff effective shear modulus 

i Index indicating position along circumference 

I moment of inertia 

I r moment of inertia of ring cross section and its effective 

width of sheet for bending in its own plane 




4 


•^strj, 

I str t 


A 

* 1 , 1=3 >*4 
L 

L 1 

m + 1 
M 


n 


Pl,P 2 

p cr 

p i 

Q 


r 
■ R 


S 

t 

U 

U r 

U 8 h 

u str r 

u str t 


NACA TN No. 1436 

moment of inertia of stringer cross Bection and its 
effective width of sheet for tending in radial 
direction (about a tangential axis) 

moment of inertia of stringer oross section and its 
effective width of sheet for tending in tangential 
direction (about a radial axis) 

index indicating position along axial direction 

trigonometric functions of cp, n, a, and t 

length of wave in axial direction 

distance between adjaoent rings 

number of ring fields involved in failure 

applied bending moment $ function of n ; a, and b appearing 
in strain energy of bending- in rings 

parameter defining wave length in circumferential 
direction 

polynomial functions of a and b 
maximum compressive force acting at buckling 
force carried by the ith stringer at buckling 
funotion of x appearing in shear strain energy 
radius of cylinder 

function of <p, n, a, and b appearing in shear strain 
energy 

number of stringer fields involved in failure 
total number of stringers in cylinder 
thickness of sheet covering 
strain energy 

bending strain energy stored in rings 
shear strain energy stored in sheet 
radial bending strain energy stored in stringers 
tangential bending strain energy stored in stringers 
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torsion strain energy stored in stringers t 
effective width, of sheet 

radial displacement of a point on a ring or a stringer 

tangential displacement of a point on a ring or a stringer 

work done "by applied forces 

axial coordinate 

angle subtended hy cutout 

coefficients used in calculation of shear strain in a 
panel due to displacements of its comers 

shear strain 

shift of neutral axis from horizontal diameter 
normal strain in a stringer 

* 

maximum compressive strain at "buckling 
angular coordinate 

THE rCEiTECHED SHAPE 


The shape of the "bulge at "buckling is determined mainly "by the 
radial deflections. The following expression was chosen to represent the 
radial deflections: 

w r = a^ sin 2 («x/L) 

= sin 2 (?ix/L)^a 0 + a^_ cos ncp + ag 008 2*19 + a 3 008 3^9 

+ "b^ sin ncp + "bg sin 2ncp + "b^ sin 3 n 9) (1) 

provided that 


0 < 9 < ( rt/n) 


(la) 


and 

(lb) 


V r = 0 
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when cp > ( jr/n) 

The notation and the sign conventions are shown in figure 2 . 

The deformations of the rings were assumed to he Inextensional. 
The condition of inextens tonality is 


v r - 



( 2 ) 


Equations (l) and ( 2 ) 


determine the tangential deflections as follows: 


w t = sin 2 {rac/L) [- a 0 <p - (aq jnj sin np - (a g y^n) sin 2ncp - ^a /3n) sin 3ncp 
+ cos nqp + (h £ y / 2n^ cos 2ncp +- cos 3ncp] (3) 


provided that 

0 < cp < (jr/n) (3a) 

Integration of the bracketed expression in the right-hand member of 
equation (l) yields an integration constant that was omitted from the 
bracketed expression in the right-hand member of equation (3) . The 
physical meaning of this constant is a rigid -body rotation of the ring. 
Moreover , ... 


w t - 0 (3b) 

when cp > (jr/n) 

If it is required that there be a smooth transition between the 
bulge and the nondistorted part of the cylinder at cp = ( ir/n) , the fol- 
lowing conditions must be satisfied: 

The tangential displacement must-vanish, that is, 

w^ = 0 ( 4 a) 


when qp = (ir/n) 

The radial displacement must vanish, that is, 

w r = 0 ( 4 b) 
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when q) = (jr/n) 

There must "be no abrupt change In the direction of the tangent, 
that is, 

&Wr/^P = 0 (4c) 

when q) = (ir/n) 

There must "be no abrupt change in the curvature, that is, 

= 0 (4d) 

when 9 = (it/n) 

The four conditions (equations (4a) to (4<£j) establish four relations 
between the Fourier coefficients and make it possible to express any four 
coefficients by means of the remaining three. If a Q , a^, and b^_ are 

retained as the basic parameters, the following four equations are 
obtained: 

&2 = (8/5)a-L - (9/5) a Q 
a 3 « (3/5) a 1 - (4/5)a Q 
b 2 - (l 6 / 5 )b 1 + ( 18 / 5 ) m Q 
b 3 = (9/5)^ + (12/5) m 0 


( 6 ) 

( 7 ) 


with the notation 

( a l/ a o) " a 

and 

cv a J = ' D 

and after substitution of equation ( 5 ) in equations (l) and ( 3 ) , 

w r = a^ sin 2 (jrs/L) 
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where r 

= 1 + a coe up + (l. 6 a - 1 . 8 ) cos 2 ncp + ( 0 . 6 a - 0 . 8 ) cos 3 nc P 

+ b sin ncp + ( 3 . 2 b + 3.6jt) sin 2 ncp + (l. 8 b + 2.4-ir ) sin 3 ncp"] (7a) 

and 

w t . = a Q kg sin^(jrx/L) ( 8 ) 

where 


k^ = (l/n) j”-ncp - a sin ncp - (l/ 2 )(l. 6 a - 1 . 8 ) sin 2 np 

- ( 1 / 3 ) ( 0 . 6 a - 0 . 8 ) sin 3 ncp + b cos ncp + (l/ 2 ) ( 3 . 2 b + 3 * 6 *) cos 2 ncp 

+ (l/3) (1.8b + 2.4 jt)cos 3 ncp] (8a) 

Equations ( 7 ) and (8) are valid, provided that 

0 < cp < (rt/n) (8b) 

When cp is greater than ( v/n) , -the deflections are assumed to vanish. 
Typical examples of the deflection patterns in the plane of the ringp 
are shown in figures 3 and 4 


CALCULATION OF STRAIN ENEROY 
Strain Energy Stored in Rings 


The strain energy stored in any one ring is 

pn/n 

TJ = (1/2) |(EI), /r3l 


jT /a [v r + ^ r /3q>| 


dcp 


(9) 


If the value of w r is substituted from equation (7) and the strain 
energy Is summed up over all the rings, the following expression Is 
obtained: 


U r = (1/2) (a^/r 3 ) ZI(El) r sin 4 ( rtXj fL) fJ [\/ a %^ 2 )] 


dcp (10) 


where m is the total number of rings included in the wave length. The 
Integration yields a result in closed form. The same is true of the 
summation if all the rings have the same bending rigidity— In such a 
case the total strain energy U r stored in all the rings becomes 

U r = (3/16) (a2/r 3 )(El) r (m + l)M (H) 
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when 


and 


■when 


m > 1 


U r = Cl/2)(a^/r3)(EI) r M 


m = 1 


(11a) 

( 12 ) 

(12a) 


where 
nM 


+ a 
+ 


■ |* + 10 . 053096^1 - 9 n 2 ) + 206.01005(1 - 4 n 2 ) 2 + 90.303387(1 - 9 n 2 ) 2 
- l 8 . 095573 (l - 4 n 2 )(l - 9 a 2 )] 

9 . 047786 e(l - 4 n 2 ) 2 - 1.5079645(1 - 9 a 2 ) 2 
30.159289(1 - n 2 ) (l - 4 n 2 ) + l 8 . 095573 (l - 4 n 2 )(l - 9 n 2 ) 

*[4(1 - n 2 ) + 2.4 (l - 9 n 2 j + ll 3 . 69784 (l - 4 n 2 ) 2 + 42.636690(1 - 9 a 2 ) 2 
2 . 4 U - n 2 )(l - 4 n 2 ) - 3.68(1 - 4 n 2 )(l - 9n 2 ) 

+ a 2 [1.5707963(1 - n 2 ) 2 + 4 . 0212386 (l - 4 n 2 ) 2 + 0.5654867(1 - 9 a 2 ) 2 ] 

+ "b 2 ^..5707963(1 - n^) 2 + l 6 . 084954 (l - 4 n 2 ) 2 + 5.08938(1 - 9a 2 } 2 ] 

+ ab j^ 6 . 4 (l - a 2 ) (l - 4 n 2 ) + 3.84(1 - 4 n 2 ) (l - 9 a 2 )_| ( 13 ) 


Strain Energy Stored in Stringers 

The strain energy shored 1 ^ the stringers "because of "bending 
in the radial direction is 


U 


str 


( 1 / 2 ) (El) 


str 


Jo 




dx 


( 14 ) 


where the summation is extended over all the stringers contained in the 
"bulge. Substitution and integration yield 


U str r =51 ( 1 / 2 )( SI )str r a o f cos 2 (2*x/L) dx 

uo 

= (^/i^a 2 5(ki 2 (EI) s 


str T 


( 15 ) 
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Both k^ and (El) g.^ vary from stringer to stringer, (ll) str 

"because the effective width of sheet to "be added to the stringer section 
changes. For this reason the summation has to "be evaluated numerically. 

The strain energy stored in the stringers "because of "bending in the 
tangential direction is 


^str-t. “ 


(1/2) (E 


dx 


( 16 ) 


where the summation is extended over all the stringers contained in the 
"bulge.. Substitution and integration yield 


a B tr t - («*‘A 3 )»12 I )s 3 (II > 


str t 


(17) 


Because "both and (El) vary from stringer to Btringer, the 

summation has to "be evaluated numerically. 


The strain energy stored in the stringer "because of torsion is 

, ,2 I"/Co \ ,12 

U t » 2_ ( 1/2)00 J ( 1/r) 


(18) 


(d 2 v r ^(Socdcp)J 

In this equation (l/r) ^8 2 w. r y(^) is the unit angle of twist of the 

stringer, and the summation is extended over all the stringers contained 
in the "bulge. In the expression for the Saint -Tenant torsional rigidity, 


C = 0.14a 


( 18 a) 


"because the test specimens were provided with square section stringers of 
edge length a. Differentiation gives 


( s \)/ ( 


(5x5cp) = a 0 k 4 (it/L) sin (2itx/L) 


(19) 


where 


= n|j- a sin ncp - (3.2a - 3-6) sin 2nsp - (1.8a - 2.4) sin 3 nC P 

+ "b cos ncp +--(6.4"b + 7>2n) cos 2ncp + (5.4b + 7 .2 it) cos 3^$] (19a) 
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Hence the strain energy of torsion is 

U t = (£/V)a| [(GC)/(lJ)] XI k 4 2 (2°) 

where the summation includes all the stringers contained in the "bulge . 
Because the variation of the torsional rigidity with effective width is 
negligibly small, the factor GC was written "before the summation sign. 
The summation was carried out numerically. 


Strain Energy of Shear Stored in Sheet 

The shear strain energy in the panel is taken as the average 
effective shear modulus G- multiplied "by the square of the average 

shear strain 7 in the panel. The value of 7 is calculated from the 
displacements of the four comers of the panel. The total strain energy 
of shear stored in the sheet then is 

' (V2>n (21) 


The effective shear modulus depends upon the geometrio and mechanical 
properties of the panel and the average strain therein. Its value was 
taken from the empirical curves established earlier at Polytechnic 
Institute of Brooklyn Aeronautical Laboratories and presented in 
figure 2k of reference 6 . 

The average angle of shear y was calculated from the equation 

7 = ( C V L l)( Wr I,J ” Vr i+l,«J “ Vr i, j+1 + Wr i+l,J+l) 

(- j + v ti + i, j ” W ti, j+i “ Vt i+1, J+l) (22) 

where the first subscript refers to the circumferential location of the 
corner of the panel and the second to the axial location, as shown in 
figure 5. The values of the factors a^, and a^. were calculated from 
the equations 

Or = ( 1 / 10 ) (d/r) = ( 1 / 10 ) (2*/3) ) 
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Substitutions yield 


Ugh = (1/2) (td/L 1 )G 0 ^ ( G eff/ G o)i R l 


(24) 


. „ number of stringer fields involved in the bulge, Q is 

^ H ” - 1 Tto sum ' 

mation Q gives a result in closed form as follows. 


m 


J=0 




^ — V 

“ jsin 2 [(*3)/(m + 1 )] - Bin 2 [*(J + D/(« + ^]j 

|(2rt)/(m + l)j 


10 (l/4) (m + l) s 


provided that 


1 - cos 


m > 1 


(25) 

(25a) 


(25b) 


(26) 


When m « 1, 

51 <4 - 2 

The meaning of the symbol R is 

B - [«r(l%i - *l,l+l) - °t(*3,l + *3,1+1)] 2 

The values of Bj. u * 3 ,i> “* *3,1+1 •» oMalned from those of *i 

and k 3 (equations (7a) and (8a)), respectively, by replacing the angle q> 

by 2*i/S or 2*(i+l)/3. As examples, k^ and- k 3,i 8X6 llstedt 
k = 1 + a cos (2«ni/s) + (1.6a - 1.8) cos (4*ni/S) 

+(0.6a - 0.8) cos (5itni/S) + b sin (2rtoi/S) 

+ ( 3 .Sb + 3.6.) Bln (W/S) + (1.8b + 2 . hit) aln (Sml/W < 2Ta ' 
t _ - (Ejmi/S)- a sin (S«ni/3) - (l/2)(1.6a - 1.8) aln (lumi/fe) 

- (1/3) (0.6a - 0.8) sin (6jmi/3) + b cos (2*ni/3) 

+(l/2) (3-2b + 3.6a) cos (W/3) + (l/3) (l*8b ♦ S.h*> COS (6«ni/S) 


the summation of the B quantities Indicated in aquation (24) ™ 
carried out numerically. 
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WORK DONE BY EUTERNAL FORCES 


As was mentioned in the INTRODUCTION, the strain was assumed to 
"be distributed linearly over the sections of the cutout portion of the 
cylinder. The force in each stringer was calculated as the product of 
the strain, the elastic modulus, and the cross-sectional area of stringer 
plus effective width of sheet. These forces were assumed to he the 
external forces applied to the stringers at the end of the axial wave 
length. Because of the reduction in the effective width of sheet on 
the compression side of the cylinder the neutral axis in bending is 
shifted toward the tension side . This shift was calculated and taken 
into aocount when the forces acting upon the ends of the stringers were 
determined. 

* 

The distance between the points of application of the forces 
shortens when the stringers bend during the buckling process . This 
shortening multiplied by the force is the work dome by the force . The 
total external work is the sum of all the work quantities calculated for 
the individual stringers: 

w = ( 1 / 2 ) [(K faT + (K/Sl (s8> 

where Pj is the external force acting upon the ith stringer and the 
summation is extended over all the stringers contained in the wave length. 
Substitutions and integration yield 

sin 2 (2*x/L) dx 

- (lA)a„(^A.)P or ^ (V^rX *!,! 2 + 

where P cr is the force acting in the most highly compressed stringer. 
The summation was carried out numerically. 

CALCULATION OF BUCKLING LOAD 


The buckling condition is 


W = (1 A)a®Z>l + *3 )£ 


^r + U Btrr + U strt + U t + U sh = W 


( 30 ) 
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where the values of the quantities must "be taken from equations (ll) or 
(12), (15) , (17 ) , (20), (24), and (29)> respectively. Equation ( 30 ) was 
solved for F or , which is a multiplying factor in the expression for 

W, and minimized by means of the following procedure: 

A value of n corresponding to a circumferential wave length 
extending over an integral- number of stringer fields was first assumed, 
and an integral number was chosen for m + 1 , the number of ring fields 
included in the axial wave length. On the basis of these tentative 
values, M, k-^, kj, and k^ were computed. Next, P Qr was assumed. 

This assumption permitted the calculation of the effective width of 
sheet and consequently of the moments of inertia of the stringers and 
made possible the determination of the values of G'&ffJ&o from the graph. 

The summations were then carried out. Substitution of the results in 
equation ( 30 ) resulted in a polynomial p-^ of the second degree in a 

and b In the left-hand member and another polynomial p 2 of the second 

degree in the right-hand member, P 2 multiplied by P cr . Solution for P 
gave the fraction 


p 1 (a,b) 
P cr “ p 2 (a,b) 


(31) 


This expression for P Qr may be minimized with respect to a and b 
by setting 

v s !l s = 

or P 2 dpg/db 

The partial differential coefficients of p^ and p 2 are linear 

functions of a and b. Equation ( 32 ) represents three connections 
between P cr , a, and b. They were solved by a rapidly converging 

trial-and -error method. First, a value was assumed for P and a 

7 C T* 

and b were determined from the two linear equations . Then the values 
of a and b vere substituted into the quadratic expression for P or . 

The procedure was repeated with the aid of new assumptions for P cr until 
the value obtained was close enough to the assumed value . 

When the value of P „ obtained in these calculations differed 

cr 

substantially from the_ value assumed at the outset, the moments of inertia 
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and the effect ire shear modulus had to he recalculated and the entire 
procedure repeated. All the calculations were carried out for a number 
of different choices of n and m. The buckling loads corresponding 
to these different values were compared, and the smallest one was con- 
sidered as the true buckling load. Details of the procedure may be seen 
from the numerical example in the appendix. 


COMPARISON OF THEORY AND EXPERIMENT 


Numerical calculations were carried out for the two arrangements 
of stringers and three circumferential sizes of cutouts investigated 
in the experiments described in reference 4. Typical buckling patterns 
obtained in the calculations are shown in figures 3 and and details 
of the .numerical results are presented in table 1. Theoretical and 
experimental bending moments at buckling are compared in figure 6 . 

The theory predicted bending moments at buckling which were, with 
one exception, consistently higher than those obtained in the experi- 
ments. Moreover, the deviations between theory and experiment increased 
systematically with decreasing circumferential length of the cutout. 
Strain-energy calculations are known to yield too high buckling loads 
when the deflected shape assumed differs from the actual shape of 
distortions . 

The circumferential wave length was predicted by theory with 
satisfactory accuracy. In the axial direction the theoretical wave 
length is greater than -the one observed. The deviation is slight in 
the case of the l 6 -strlnger specimens and large In the case of the 
8 -stringer specimens. Small changes in the axial wave length, how- 
ever, have little effect upon the buckling load. 


CONCLUSIONS 


A strain-energy theory has been developed for the calculation of 
the buckling load In general Instability of circular reinforced mono- 
coque cylinders which have a side cutout and are subjected to pure 
bending. When the theory was applied to the test cylinders of part IV 
of the present series of Investigations, the following percentage 
deviations from the experimental values were obtained: 54.8, 32 .4, 

and l 6.1 percent for the 45°, 90 °, and 135 ° cutouts, respectively, of 
the 8 -stringer series; and 47*4, 30 . 8 , and — 3.2 percent for the 45° # 
90 °, and 135 ° cutouts, respectively, of the l 6 - stringer series. " 



1 6 


NACA TN No. 1%36 


The authors "believe it would "be desirable to continue the theoretical 
investigations on the basis of more refined assumptions for the deflected 
shape at buckling. Strain -energy calculations are known to yield too 
high buckling loads when the deflected shape assumed differs from the 
actual shape of distortions. More experimental work is also needed for 
a better understanding of the general instability phenomenon of reinforced 
monocoque cylinders having a side cutout v 


Polytechnic Institute of Brooklyn 

Brooklyn, N. Y., December 26, 19^6 
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appendix 


HUMERICAL EXAMPIF 

Details of the numerical -work performed in connection with the 
determination of the "buckling load for a test specimen are shown. The 
cylinder considered is cylinder 39 of reference 4 . Pertinent data to "be 
referred to are listed as follows: 


Radius , r, in. . . ............ 

Distance "between adjacent rings , L]_, in 

Width of panel measured along circumference, d, in.. . . 

Angle subtended by cutout, a, deg ...... 

Total number of stringers in cylinder, S 

Stringer cross section, in.. ...... 

Ring cross section, in - 

Young's modulus, E, psi . 

Shear modulus, G, psi ......... 

Shear modulus of sheet covering at zero compressive 
• load, G Q , psi 

Thickness of sheet covering, t, in 

Moment of inertia of ring cross section and its effective 
width of sheet for bending in its. own plane, I r 
(( 1 / 2 ) ( 3 / 8 ) [( 1 / 8 ) + (0.01213), 


. 10 
6.1+29 

7-854 
. ' 45 


8 


3/8 x 3/8 
3/8 x 3/8 


10 .5 x 10 6 
3.9 X 10 6 


3-9 x 10 6 
. . 0.012 


80.35 x lo- 6 


Once a value of s, the number of stringer fields involved in 
failure, is chosen, it is possible to reduce the expressions denoted by 
M, k^, k^, and k^ defined by equations (13) , ( 7 a), (8a), and ( 19 a), 

respectively, to arithmetic polynomials in a and b. The one for M 
is obtained by simply inserting the values of n = (s/2s) and its powers 
into equation (13) . For s = 3, the result is: 

MX io -5 = 0.0020877a 2 + 0.0131007b 2 + 0. 19654065 


+ 0.002868l5ab + 0.0084369a + 0.1010576b 


( 33 ) 


For the evaluation of k^, k^, and k^ it is convenient to set up a 

tabular arrangement. For s = 3 this arrangement takes the following 
form: 



CO 


Jt 

JL 

Constant 

oaa (iti/j) 

cos ( 2 id/ 3 ) 

cos ( 3 * 1 / 3 ) 

2 * 1/8 

sin (jd/ 3 ) 

sin ( 2 jd/ 3 ) 

sin ( 3 id/ 3 ) 

0 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 

.5 

-5 

-1 

.785398 

.8660255 

.8660255 

0 

2 

1 

-.5 

-•5 

1 

1.570796 

.8660255 

-.8660255 

0 

Multi- 
pliers 
for k^ 

1 

a 

1 . 6 a- 

1.8 

0 0 

® jr 

■ 


b 

3 . 2 b+ 

11.3097336 

1 . 8 b+ 

7.5398224 

Multi- 

pliers 

forkj 


0.75b 

1.2b+ 

4.2411501 

0.45b+ 

1.8849556 

-1 

-0.75a 

- 0 . 6 a+ 

O .675 

- 0 . 15 a+ 
0.2 

Multi- 
pliers 
for k^ 


1.333 •• *4 

8.5333 ...D+ 
30.1592894 

7.2h+ 

30.1592894 


-1.333. -.a 

-4.2666. . .a+ 
4.8 

-2.4a+ 

3.2 


For tlie value of the index 1 = 0 denoting the position of the stringer at the edge of the cutout, 
the polynomial for ^ in a and h is found, "by multiplying the expressions appearing in the first 

row "below the double line (labeled multipliers for k^) by the numbers in the corresponding column 

listed in the first row (labeled 0) and adding like quantities of the results of the products. In a 
similar manner the polynomials for i = 1 and 1 = 2 for are obtained, as well as the three 

values each for k^ and k^. The results are presented in tabular form as follows: 
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TABLE 2 


■ 

k i 

k 3 

\ 

b 

3.2a-1.6 

2.4b + 6.1261056 

17. 0666... b +60.3185787 

B 

- 0.9a + 3.673071b + 12.4945177 

- 1.1691344a - 0.675b - 4.2063614 

- 4.849728a - 10.8b - 41.0820117 

E 

- 0.7a - 1.9052561b - 8.6945177 

- 0.1299038a - 0.525b - 2.3909827 

2.5403415a + 2.2666.. b + 10.9227221 


The functions (k-^ i - k-^ 1+ i) and (kg i+1 + kg can he determined with the aid of table 2 by 

simply subtracting or adding the polynomials of adjacent rows in the first two columns of the table. 

It must be remembered that k^ = 0 and kg =» 0 when i - 3 since the deflections have been assumed to 

vanish at the third stringer, where i * 3. If is taken as O.O785 and as -0*5, the function 

E = (k, . - k. -i 1 1 )a - (k„ . ,, + k . )a, becomes for each field: 
x,x-ro. r 3,1+x 3,1 t 


TABLE 3 


ith field 

E 

1 

-0.2625532a + 0.576826b - 0.147112 

2 

-0,6652272a - 0. 164686b - 1.6344854 

3 

-o.ll99299a - 0.4121398b - 1.878358 


H 

VO 


TN No. 1436 

















Before the stringer bending strain energy, shear strain energy, and external work can he 
evaluated. it is necessary to assume first a value of the buckling load so that the numerical 
values of the stringer moments of inertia, effective shear modulus , and effective area may he 
determined. Since the last three quantities are functions of the normal strain acting in the 
axial fibers at buckling, a value of the critical strain is assumed. Also in order to locate 
the position of the neutral axis at failure, the critical strain must he known. In the present 

example, the critical strain is guessed to he 21 X 10 -4 . Then the shift of the neutral axis, 
calculated by taking fiyst mcsnentB of area, is found to he 0.07, expressed in percentage of 
the radiuB. The following table contains the af ore-mentioned items for 6 » 21 x 10 “4 

and 6/r = 0.07 ““ 


TABLE 1* 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

D 

€ 

2v 

^str r 

I str t 

6 

3.3 


A eff 


0 

8.884 x lo -4 

4.290 

22 .15 X 10 -i t 

435 X 10"4 

2.69 X 10 4 

0.633 

c.1664 

0 .0704 

1 

19.506 

2.648 

23.85 

200 

5.91 

.419 

.1724 

.16014 

2 

19.506 

2.648 

23.85 

200 

5-91 

.419 

.1724 

.16014 


Column (l) refers to the stringer station. Column (2) , the strain at these locations, is 
directly proportional to the distance from the neutral axis, since a linear strain distribution 
is assumed. Column (3) is the effective width of curved Bheet calculated from equation (30) 
of reference 7* Columns (4) and {5) give the moments of inertia of the stringers plus their 
effective width of curved sheet. These are calculated from equations (34) and (36) of 
reference 7, with special considerations for the edge stringer because it has effective width 


1 

i 
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on only one side. Column (6) indicates the ratio of the actual strain 
in the sheet when the monocoque cylinder "buckles to the "buckling strain 

of a panel of sheet of 3*30 X 10 “^" as given in reference . 6. These values 
are needed to obtain the percentage reduction in the value of the shear 
modulus recorded in figure 2 k of reference 6 and presented in column (7) • 
Adding to the cross-sectional area of the stringer, 0 .lk 0625 square inch, 
the area of the effective width of sheet which is 0.012 times a value of 
column (3) yields a value of effective area shown in column (8). The 
entries in column (9) can "be computed with the aid of columns (8) and 
( 3 )- 


The next step in the calculations is to assume a value of m + 1 
the number of ring fields involved in failure. It is then possible to 
write the equation for the buckling condition in terms of the parameters 
a and b . It is convenient to multiply the numerator and denominator 

terms of equation ( 31 ) by £(m +l) /e] X lcA and to solve for X 10^ 

instead of this procedure requires that each strain energy be 

multiplied by £(m+l) /efj X 10 ^ and that the external work be multiplied 

by ftm+l) /eTJ X 10 \ When m +1 = 11 , from equations (ll), ( 15 ) , ( 17 ), 

(20) , (2k) , and (29) , there results : 

U r pm +1 )/eJ X 10 1 *" = a^ [(3/16) (11) 2 (80. 35 X 10 _6 )/10 3 ]m X 10^ 

2 

= a Q 1822.9k (MX 10 - 5 ) (3k) 

where M X 10 - -'’ is given in equation (33) . 

U a tr[“ +D /4 X 10 1 * = a®[ K 4 /fe.lie9) 3 (ll) a ][lI + 5Z lc f I e tr t ] * ^ 

“ 4 11 [( 30-298 + ( 30.298 ( 35 ) 

where the indicated summation can be evaluated by taking the Bum of the 
products of the squares of the polynomials for both k-j_ and appearing 

in table 2 and 3O.298 times the entries for corresponding values of i 
listed in columns (k) and (5) of table k. For this purpose it is help- 
ful first to calculate and record the quadratic polynomials in a and b 



which represent the squares of k, and k, for each value of i . These expressions are not shewn 


here . 


4 


(m + l)/E 


h 2 

X ID = 8l 

0 


= a 


10) 2 ( 6 . 1*29) 

(0.039^7)^ XL kf^ 


( 3 . 9 /IO. 5 ) (0.002769) 


(?<) 


X 10 


(36) 


■where the summation is performed by squaring and adding the values of k^ given in table 2 for each 
value of i . 

I 

n»h[t* + x - % [(0 .002) (t - 054) / ( B) (6 -itag)] < 3 - 9/io . 5) Kn) 2 (i - cm fjj 


^( G e«A)^] * lol> ‘ V 


(37) 


iS 


where the summation can he deduced by squaring each polynomial of table 3 , multiplying the result by 
13°*57 times the proper value of (& e ff to be found in column ( 7 ) of table 4, and adding all such 
products. A table containing the terms in is needed. It is not given here. And finally, 

Vf [(« | 1 ) £] x lO 1 * - a* /(i, (6 ,te 9 )] Vt x 1 <A h_ Vr («/V«t)(l i + 4 j 

= a o'mi Xwlt |^ |(0 ■3 8 379)A e:ff ( , /«„ ix )](k| + k|)| ( 38 ) 

where the summation is carried out by finding for each value of i ^he sum of the squares of k-j_ and 

k, . These quantities are evaluated in connection with the determination of the stringer bending strain 
energy, multiplying this sum. by 0. 38379 times the corresponding value appearing in column ( 9 ) of 
table k, and adding all such products. 
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The value of € X 1CT is obtained, by adding equations (3^) , (35)? (3^) , and ( 37 ) and dividing 
max 

this sum by equation ( 38 ) divided by ^ x 10 4 . If the foregoing procedure is followed, the numer- 
ator of the ratio of the two quadratic polynomials in a and b to which ^ X 10^ is equal can 
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In the present numerical example e raea: x 10^ has been assumed to 

"be 21 for the calculation of the position of the neutral axis, effective 
widths of -sheet, moments of -inertia, and effective areas. The same value 
is substituted into equations (40) and (4l) . The values of a and h are 
found to "be 


a = - 3 A 349 
b = -1.74254 


(42) 


When these results are inserted into equation (39)j ^ X 10^ becomes 
22.44. Experience has shown that the difference between the value 
6 max = 22.44 X 10-* obtained and the value c ma;x; = 21 X 10 -i+ assumed 


is small enough to make a repetition of the calculations unnecessary. 
The changes In the values of the effective widths, moments of inertia, 
location of neutral axis, and effective- shear modulus resulting from 

the increase In from 21 X 10“ ^ to 22.44 X 10" ^ would have a 

negligible effect upon the calculated buckling strain. Hence 

22.44 X 10 -J+ may be taken as the critical buckling strain. 


The critical moment is computed from the relation 

- n pa 


where P Is the force In a stringer at buckling, d Is the distance 
of the stringer from the neutral axis, and the summation must be taken 
over the entire cylinder. Experessed in termB of strain, this relation 
is 


M, 


cr 




A eff d 2 


(43) 


since 


e 


e max d / (r + S > 


Eor a critical strain of 22.44 X 10“\ M cr becomes 174,960 inch-pounds 
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s 

a 

(deg) 

Mtheor 
(in. -lb) 

^exp 
(in. -lb) 

Percent 

Difference 

8 

45 

90 

135 

174,960 

142,890 

116,070 

113,000 

107,900 

100,000 

54.8 

32.4 

16.1 







1 6 


45 

90 

135 


342,680 232,500 
225,690 172,6oo 
164,470 169,900 


47.4 

30.8 

-3.2 
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SECTION A- A 
NUMBERS IN PARENTHESES 

REFER TO FIELDS 


Figure 2.- Deflected shape. 
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n = 
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n 

= 1* 

a = 

- 2.7739 
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= -2.3523 

b = 

- 2.7464 

b 

“ - 3.1193 
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Figure 5.- Deformations of panel corners with notation for shear strain calculations. 



